


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1978 


The simulation of gas turbines by a state of 
the art analog device. 


Shepherd, William M. 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/18546 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


THE SIMULATION OF GAS TURBINES BY A STATE 
OF THE ART ANALOG DEVICE 


William M. Sheperd 


7 MAR 1979 








eer ra a ee 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE READ INSTRUCTIONS 


BEZCRE COMPLETING FORM 
2. GOVT ACCESSION NO] 3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (and Subisitie) $. TYPE OF REPORT & PERIOD COVERED 


eres SIMULATION OF GAS TURBINES BY A STATE THESIS 
OF THE ARI ANALOG DEVICE 6. PERFOAMING ORG. REPORT HUMBER 


7. AUTHOR(ae) - CONTRACT OR GRANT NUMBER( 


SHEPHERD, WILLIAM M. 


9 PERFORMING ORGANIZATION NAME AND AODRESS - PROGRAM ELEMENT, PROJECT, TASK 
AREA &@ WORK UNIT NUMBERS 


MASS. INST. OF TECHNOLOGY 


CONTROLLING OF FICE NAME AND AODORESS 12. REPORT DATE 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA, 93940 ES 


MONITORING AGENCY NAME & ADORESSE(// different frem Controlling Office) 


DECL ASSIFICATION/ DOWNGRADING 
SCH EOULE 


DISTRIBUTION STATEMENT (ef this Repert) 


APPROVED FOR PUBLIC RELEASE; DISTRIBUTION UNLIMITED 


OISTRIBUTION STATEMENT (of the ebetract entered in Bleck 20, If different frem Report) 





16 SUPPLEMENTARY NOTES 


KEY WORDS (Continue on reveree eide if neceeeary and identify by biock number) 


GAS TURBINES; ANALOG DEVICE; SIMULATION OF GAS TURBINES 


ABSTRACT (Continue on reveree aide If necessary and identify by bleek mumber) 


SEE REVERSE. 





OD an 7s 1473 eolvion oF | Nov e815 OmsoLeTe UNCLASS 
(Page 1) SECURITY CLASSIFICATION OF THIS PAGE (Phe Dave Snive 





ay ea 


THE SIMULATION OF GAS TURBINES BY A 
STATE OF THE ART ANALOG DEVICE 


by 


William McMichael Shepherd 


Submitted to the Department of Ocean Engineering on 
12 May 1978 in partial fulfillment of the requirements for 
the Degree of Ocean Engineer and the Degree of Master of 
Science in Mechanical Engineering. 


ABSTRACT 


A simulator was designéd to model the operation of a 
Single shaft gas turbine engine. The engine characteristics 
were represented by a third order non linear mathematical 
model implemented with analog computation. The use of 
"state of the art” integrated circuitry allowed for a 
considerable reduction in the space and power required for 
this device compared with conventional analog methods. 
A hard-wired desk-top machine was fabricated to simulate 
the real time dynamic behavior of a gas turbine prime 
mover. 


Thesis Supervisor: Henry M. Paynter 
Title: Professor of Mechanical Engineering 


Thesis Reader: A. Douglas Carmichael 
Title: Professor of Power Engineering 


7195835 








THE SIMULATION OF GAS TURBINES BY A 
STATE OF THE ART ANALOG DEVICE 


by 
WILLIAM MCMICHAEL SHEPHERD 


B.S.A.E., U.S. Naval Academy 
(1971) 


SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
DEGREE OF 
OCEAN ENGINEER 
AND THE DEGREE OF 
MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


AT THE 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
May 1978 


(C) William McMichael Shepherd 





Ee 


THE SIMULATION OF GAS TURBINES BY A 
STATE OF THE ART ANALOG DEVICE 


by 


William McMichael Shepherd 


Submitted to the Department of Ocean Engineering on 
12 May 1978 in partial fulfillment of the requirements for 
the Degree of Ocean Engineer and the Degree of Master of 
Science in Mechanical Engineering. 


ABSTRACT 


A simulator WaS designéd to model the operation of a 
Single shaft gas turbine engine. The engine characteristics 
were represented by a third order non linear mathematical 
model implemented with analog computation. The use of 
"state of the art" integrated circuitry allowed for a 
considerable reduction in the space and power required for 
this device compared with conventional analog methods. 
A hard-wired desk-top machine was fabricated to simulate 
the real time dynamic behavior of a gas turbine prime 
mover. 


Thesis Supervisor: Henry M. Paynter 
Title: Professor of Mechanical Engineering 


Thesis Reader: A. Douglas Carmichael 
Title: Professor of Power Engineering 





ACKNOWLEDGEMENTS 


The author wishes to acknowledge Professors 
H. M. Paynter and A. D. Carmichael for their assistance, 
advice, and guidance in the completion of this paper. 
Special thanks is also extended to Mr. Ken Busick of 
the Beede Instrument Company, and Mr. John Carter of 
the Henschel Corporation for providing much needed 


hardware. 





TABLE OF CONTENTS 


Page 
MRRIIG RAGS «sec cee cee tec ce dee AON 5 ee 1 
eee aie CMe. tee et NOT ee la eS cc wc et ee eee ee rates 2 
ROW ICAGGMENtS smemeetS . . Shee ss ccc cc cc ct cece ee eee eee 3 
aS Le Cie Tee eee W eee On eae 4 
i. BME OCAWEEHOMm reise e . sacs « @ juin 66 7 
A. Background ....... ee Se eee 7 
Baw ) DeSicieCameept ie. .ssc cess 6. ss cbse tee ree 9 
el, MIE GI EOC pes neice lee Ly 
iN 9 ES Gygsl J sioaeys le) 11 
B. Steady-State and Dimensional Analysis .... ll 
C. Frequency Analysis ........... ee a ee 16 
D. Dynamic Equations ............. ee arate 16 
Iho INE SUNT Oa Ke) (Cee ee 16 
PRO me EMGUSO UG) «co cic ct sce ses eer eoeces 17 
Bee MUI eGapaettanCe 6. <2. 5 ce sce ee ccc ce ee 17 
Hee State mquabions ........... eee ee 21 
F. Engine Data, Assumptions, and 
MPR ACHOTEUSEVOCS cic ccd sees veer csccencecees 24 
Si Pips IST ON | aisle cue cic cs 6 eueeie + coc cc cece 24 
2. Design Point Sepcifications .......... 25 


3. Compressor and Turbine Characteristics @25 





0 Ga 


IV. 


VI. 


ie "OPEraclOMma lL AMPUNE LICL ices ee ce we ecw eces 
Pe 46INGROUBCEION (oN ee ss sss ss Ser Oe eC ee 
B. Basic Operation ..... 5 6 Ee 
C. Inverting OP AMP ...... 5) OO eee 
D. iInverting Summer ..... WOME fees Pon thoy en cei" ie 6. 6:6 sce 
ye eT) COC ee @ laieicicls «1c cle so 6s 0 sss 00 3 6 giaietie! = ices «)s' ¢ 
Re LOG aa teminmewGeneratOYr 6. .6ss csc sew swe ewes 
Cem ATIC LOGMGENCTALOL 26 sess 6 st sec e cece eee aoe 
Raa lOg s@ommmeadclON ...<.se cece ce se seco Cen ene.e s 
Perm tedgerelie @ LON 466 6.6 0 0s ¢ sieiee 0 e Oe ne eae 
Bue tiem bractical Log Circult ....¢s.6.<ss Soe 
Ema iiclian mC LGM (seis «ac icc ole ce sess cicceseees 
Pee POGUCELON sls c 0:5 0 «s/s MeN reM Crete eile e ehec RUS ae 6) 
Be meEleCteronics ..... Peeters. ai. . ee 

ILA LGXe On) elet tl 

2. DYSYS Simulation and Scaling ..... agave 
EEO NIG Teak US Ue sic co clsele soc e se cscs cece ccceee a 
Implementation and Checkout ............ secmeaesrs 
Pere Ae CHEUIR@ Me (slices <5 4 c/s 0 cle ssc ¢ os 006000 0006s 
DPIG@CMGCTROGACLON 12.5. c cscs sce ccs secseecoces 
SPEPECCOCINOCNELAELON 26.00 c6c ec cece wcscsecvecesecs 


DeGeMeGoalwCUrCULt DESIGN ...e.cese cess ces eee 


29 
29 
30 
Syl 
32 
33 
35 
36 
36 
36 
4l 
41 
41 
41 
42 
45 
47 
47 
47 
47 


48 








Page 

vid. Conclusions and Recommendations ............. a1 
Peeetecacibatsity of thesPesigm .......2ss..... ol 

ci WU pe¥eniei) etoile tes! 0) o1 
REGEYTENCES .c.ercseses EN 00 0 OE 53 
Appendix I - Machine Characteristics .........ceeee. 25 
Appendix II - State Equations ..... ONO R666 OOO Toe 62 
Pependix WU Tcce Datal........0....-.-s0.-e-00- 66 
Appendix IV - DYSYS Simulation ............cceceeeee 93 
Appendix V - Design Drawings and Sketches .......... 97 
DZ 


ABBeNGLX Vill opeesseteatdONo G2. 5.. ose et ce tt ee 





I. INTRODUCTION 


A. BACKGROUND 
This paper describes the concept and construction 
of a hard wired analog device, pictured below, which 


Simulates the operation of a gas turbine engine. 





Recent work on gas turbine simulation has focused 
largely on digital or hybrid computing. To accommodate 
the nonlinear elements of a gas turbine simulation in this 
fashion requires a perturbation or piecewise linear 


technique. A pure analog approach offers a continuous 
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solution to these nonlinear characteristics, which facili- 
tates a Simulation in real time. 

At this time, analog computation of complex non- 
linear systems is not feasible on devices such as the EAI 
680 computer currently installed in the M.I.T. Joint 
Computer Facility. With the miniaturization of both digital 
and linear functions into the now-common integrated circuit 
or IC chip, it becomes possible to design extremely compact 
analog devices. The IC component cost of such a simulation 
is rapidly approaching the cost of purely passive electronic 
elements. 

The modeling approach to the analog simulation, 
physical modeling, follows the work of Markunas- closely. 
Physical modeling attempts to explain the engine char- 
acteristics within a framework of physical princples and 
laws which describe the functions occurring in the real 
machine. Simplicity is generally a strong point of the 
model. 

The model which has been implemented in the analog 
simulation represents a general machine "invented" by 
Markunas, who estimated the necessary characteristic 
relations between the independent variables. Markunas 
conducted digital and hybrid simulations of his "general" 
gas turbine, but a full analog implementation was not 


investigated. 
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This paper is essentially a design manual which 
traces the development of a hard-wired analog simulator. 
It is intended to give the reader a brief overview of the 
analog devices currently available, and their application 
to a computational circuit. To facilitate future work 
with the simulator, a detailed user's manual has been 
written. 

This document contains explicit details relating 
to the construction and operation of the simulator. The 
user's manual, the simulator, and associated hardware 
are under the supervisor of Professor Henry M. Paynter, 


Department of Mechanical Engineering. 


ee DESIGN CONCEPT 
The following characteristics were desired in the 
hard-wired simulator: 
l. Real time, non-linear third order simulation 
of the dynamics of a single shaft gas 
turbine engine, from start-up to overload 
conditions. 
2. Cost to be under $600. 
3. Size of the device such as to be readily 
portable and convenient to set-up and operate. 
4. Layout and operation of the simulator 


to be comparable with typical industrial 
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and military engine control consoles. 
Accuracy of the steady-state simulation to 
be considered of secondary importance 
compared with the representation of large 
scale transient effects. 

The ability to investigate various fuel 
control schemes in a closed loop control, 
and the ability to operate with different 
torque-speed load characteristics. 
Utilization of the general purpose, low 
cost IC operational amplifier in a hard 
wired configuration to accomplish all 


Mathematical operations. 








ale: 
If. MODELING 


A. GENERAL APPROACH 

The gas turbine model selected by Markunas will be 
briefly explained. The modeling approach is broken down 
into: steady state and dimensional considerations, 
frequency analysis, and dynamic equations. 

Markunas' estimated performance parameters and 
the dynamic equations are combined in state-variable 
format. These equations appear in Appendix II and form 


the basis for the analog signal flow used in the deSign. 


B. STEADY STATE AND DIMENSIONAL ANALYSIS 

Horlock> establishes non-dimensional functions in 
describing the performance of compressors and turbines. 
A set of relations are developed to describe the 
performance of a turbomachine in its simplest terms. 

In the subsequent development, gas turbine state 
points are used, numbered in the following convention, with 
total values indicated by a zero in front of the state 


point: 





Sar 


COM PRS ar. TURBINE 


INLET EXKKAQSTT 
om 
aa } 


"2 


2 s at ae 


The following symbols are used for dimensional 


analysis in the compressor: 


S ol 
Sz 
COZ 
ros 
a2 


g A 


Compressor 


Inlet stagnation pressure 

Inlet stagnation density 

Outlet stagnation pressure 

Inlet stagnation temperature 
Efficiency 

A characteristic linear dimension 


Stagnation temperature rise through 
compressor 


Rotational speed 


Mass flow rate of gas 


Dimensions 


m/_Tt < 
40) ee 
e/a 








Ps lac 


Symbol Compressor Dimensions 
vi ei Molecular wt of gas aa 
A Absolute viscosity of gas m4 /LT 
CP Specific heat at constant pressure Le 
27 Ratio of specific heats 
rR Gas constant oo 1S 
Aez Speed of sound of entering gas _/ae 
at Torque input mMia/T? 


In the compressor, if Foe Made “ee and ANo are chosen as 
variables of interest, one can write 

Fes,4, Ti, Alo = F(z, Tee ND, m Coz, & 4, Aoz ) 
Assuming a perfect gas, fixing the gas constant,f2, choosing 
“S, and Taz thus constrains “402. Ose can be expressed in 
terms of Kee loz ; the original relation becomes | 
Pin, Te, ATR = © (a, Toe NDR Kw) 
Taking non-dimensional groups, 


— Sah AVo _ ¢ ND. mR Be ta2D y8 
Sz 1 Bz2D* , Tez Alaz R207 -VTee 


With constant values, and xe can be eliminated: 





Bot Ges: "Te 
Che Seis. a CN [az 


Multiplying the last two non-dimensional parameters to 


create a new parameter, 








Ae 


Fou ere 4 NS = NID mViez IN 


cate eens OOO 
we yaa! 


Ry. Rep te oz. ic lez Be or AAD 
10 ol be that the third paramet oe 
can seen a p ye isa 


"flow Reynolds number". 
Assuming that the compressor operates above some 
critical Reynolds number, and that the dimensions of a 


particular machine are fixed, the equations reduce to: 





As yt AT Le (HN. mh Tez 
Paz ees ea NTs. Foz 


In order to present meaningful performance relations, 
compressor data should be corrected following the above 
parameters. 

Identically for a turbine,” it can be shown that 
the same relations hold and that the parameters of 
interest in turbine performance are m,n V2 Ale 
All inlet conditions assume values at state point 4, and 
all exit values are at state point 5. Torque across the 


: . ote 
turbine is !2 . As shown before, 


NTs 7 ATs 


7 
Peg ew Tha ies a 


The characteristics for a given Ren onacHine can 
now be specified as functions of dimensionless or pseudo- 


dimensionless groups: 





Compressor: 


volTez, 


oz 


faz 
Tez 


et 


on pa 


Yo2 


Turbine: 


MN Tot 
Pod 


Tas 
TO4+ 


=o 


fs NN 

= + Caz. , VEE S 
('S . 

i" ( ee co 





eae aheeeeneell 


= (For = 





rs WN. 5 
7 Pea ) (Tod 
ae 7 Sy 
ete \ Tes 








-l16- 


c FREQUENCY ANALYSIS 
Modeling of the processes which take place in a gas 
turbine can generally be classified as mechanical, fluid, 
or thermal in nature. To reduce the scope of the problem 
to a more manageable field, only the most significant 
energy domains and dynamics were considered. Markunas 
defines three dynamic regimes for the gas turbine: 
1. Low frequency, (f<—~ .1hz) where rotary 
inertias and thermal capacitances dominate 
the dynamics. 
2. Medium frequency (f<~ 1 hz) where rotary 
inertias and fluid capacitances dominate. 
3. High frequency (f 100 hz) where the fluid 
dynamics are dominant. 
Markunas established that rotor dynamics and fluid 
capacitances were the primary influences on gas turbine 
dynamics. Thermal capacitances and fluid inertias were 


neglected. 


D. DYNAMIC EQUATIONS 

1. Assumptions 

The gas turbine model is assumed to operate with a 
perfect gas in one-dimensional flow. Components are 


treated in lumped parameter fashion. 





2. 
For 


following: 
© 
Ne 
where 
Nic 
Le 


a 


3. 


sy 


Rotor Inertia 


the rotor, a summation of torques give the 


an 
—— 


angular acceleration, rpm/sec 


rotary inertia, ft-1bf-sec* 


load, compressor, or turbine torque, ft-lbf 


Fluid Capacitance 


For the one-dimensional lumped parameter fluid, the 


following symbology is used: 


fluid density, lbm/ft> 

axlal distance along duct, ft 

mean axial velocity of fluid, ft/sec 

ratio of specific heats = Ce/Cy 
specific heat at constant pressure 
specific heat at constant volume 

Mach number = ‘NX / OO __ 

Local sonic velocity = Ey SRT 
internal energy of the fluid in ft-lbf/lbm 
enthalpy of the fluid in ft-lbf/lbm 


lumped fluid volume 





=n = 


A general mass conservation equation for the one 
dimensional fluid stream may be written: 
=m Ce) > Vx = Gy 
—_ 
“Ot tke 
Integrating over a constant area duct, with fluid pro- 


perties constant at a given axial distance or dimension: 


a (fA) , BR CAN) = 


“Oo 
DA. 


ee 
Integrating over the length of the duct 


2S ¢Adx + { BCAv« dx = 0 


on 


by definition. 





a (Ee -vor ro ANS 1S 
aCe CAN | | 
§ eds is approximated by an average density, 


——— 


¢ e» Lumped fluid volume (Ver) since 


ate) = .... CW ~ ae ) 
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where the subscipts 1 and 2 identify inlet and outlet 


conditions across the engine, respectively. 


Writing a similar equation for the conservation of 


momentum neglecting viscous shear and body forces, 
CB) + RCL +P) =8 
a 


—/ over area, 


3. CRRA) +> 5, (0A A + PAY _ 


3. (CBA) + BCBS) a > 


- = es 
BME) + cee A2e + sa @) = 


Integrating over length of the duct 


no 
Be (CMA) bx rm 


“Zz “e 


+ PA 
rs ae 


~—2 a” 
+ = ‘a; ls ~ PA | ce 
Sr mL a 








4 (Bt) - PA + tue : 
ae 

Sy. 3h 
OH) = LPs eye] 








= j= 
a 3 


y = ~ 
v2 vy; = x 
Mos ae = cert 


a Es 
ra A =e 
£R)- PAP we), 
Essen 
Since total pressure, PA e(L Ne s-= at he 


Rss | Ble <PV «™M~ » 


ae 7 C+ M7)” 5h 


Assuming Y = 1.4, and M Z .5, the momentum equation can 


be written: 
tC) = +. . - Sz) 


Similarly for energy conservation with no external 


work or heat addition, 


a fe(orye >| * Be (h V4 )| ae 


“2AC 


Integrating over area and length, 


— Fea(uyst S >| eindh st) - O 


3S ZAC 2s" 





2 


ake 





= 


(vor) $e C0 a > | = AN \N\0 


where = 
Wo = \n a Nes 
—w yO 
Since a _ ~ 
Vx = ~RTM 
wesc 2 
Substituting, 


Ba ia” 2) GS CL + 3G ‘ 


T= 





de ( ae Ts | i+ SCL 2 SE (inex -treTee ) 
Y et we 
= 


where Ta = T+ (t+ BS m=) 


Cae 
With M = ay the result can be simplified to: 


J (8) =X Cartas - te Ke 
a Hs Vou 
E. STATE EQUATIONS 


The three conservation equations are now: 


1. Mass 


ee) tM — 2) 





m9 


2. Momentum 
=\ ack | > 
3 Cn) = 2S (Ri - Sz 
3. Energy 


2 (eT) as S (Auk -M-s ==) 


Using the definition of a polytropic process, 


a) 
ry = CONST, it follows that 


fo. = CCONST ) £0 
ere a 
— . Cconst)% “= Const) 6 


IS 


In order to get the mass and energy equations into more 


convenient form, 


gC) ~ Ceonst)(o) dL C25) 


28) Piaie ZT )4 rei CB) 


al 4 Ses ': Ve + a cd alk 
a € St £% as (n- ~i\conss $s 


 (& (4.7) 2. EY 1) 


y. WHECE a iS ALSUMEO TO AK srmenarion 
UI Aen ss 





= 


With the result that the equations now be more: 


CP CLAS TID 


4. (RQ) = SEE (mi - che) 
Se (™) 


Seto) = SGERR.. re - Rater) 





<A (Re fz) 
a 


Where 1 and 2 represent inlet and outlet across the 
duct. Because of lumped parameter considerations, machine 


model parameters can be substituted as follows: 


One Dimensional Machine Model 
WW, Ws 
We W4 
= Taq 
To. Tos 
Teo Vos 
Po Tox 


Substituting notation for the first and third 


differential equations, 


Fos = (OBES) [a wha] 








a) Ales 


(tyes) Bites ) “slog +i -¢ 
oe = a Foe VoL iis + reba) “4 


Introducing the torque equation as the third system 


equation, 


where 

hay = fuel lower heating value, Btu/lbm 
> 

We = fuel flow, lbm/sec 

ke 

\ az, = load torque, ft-lbf 


rez, Tad, AnNO Ne are chosen as variables of 
convenience. 
For an isentropic process, n = 1.4. Markunas felt 
that n-l/n should be replaced by 1.0 in the temperature 


equation. 


F. ENGINE DATA, ASSUMPTIONS, AND CHARACTERISTICS 

i 2acroauction 

In Markunas' analysis a single shaft machine was 
invented from a number of sources. Several real world 


Machines exist which have design parameters close to 





rs ke 


Markunas's engine. The recent Garrett GTPF 990 engine in 


its locked shaft mode is a reasonable example. 


2. Design Point Specifications 


Pressure Ratio SiO) 
Compressor Isentropic Efficiency 1280°F 
Turbine Isentropic Efficiency 88% 
Compressor Discharge Mass Flow 100 lbm/sec 
Inlet Pressure 14.175 psia 
Inlet Temperature 80°F 
Rotor Speed 7200 rpm 
Output 6948 HP 


Markunas plotted characteristic curves for his 
compressor and turbine, which appear in Appendix I. 
Markunas then mathematically "fit" each characteristic 
curve. These approximations follow the functional 
relations established for non-dimensional and pseudo 


dimensionless groups. 


3. Compressor and Turbine Characteristics 


Markunas defined the following variables and char- 


acteristics for his compressor: 


PRO = z/Rz pressure ratio 





SoG 


TRO = Tox / loz temperature ratio 





= Wa Tee yz, Voz scaled compressor 
=" as he mass flow 
‘DESIGN 
KNC = Noe (Ss p) scaled compressor 
= — speed 
loz STZ Lice 
xMC = ae (rn3NTez ) 
\Tez Pes2 OFSiIGN 
, <= 
AC = .75 CXNS) + 27 CKNC) 


as 


ic qi 
Sc = SO C%NIC ) + 3.0 CXNC ) 


s.8 
Sen = sis (2 — XNC) 
a 
Cc = (wi . gy - OXNC —9) 


4-4 
EC = L.o + XNC + 3.6 (KNC) 


eC = 4.0 \e PRC £ EC 


eG iF PRO weEC 





ay de 


MWC = AC] L — (rmesS) 


EC 

EFFC = Oct : ae 
aoe 

TRE = |r (recy 


EFFC 


Similarly, for the turbine 


rev a Yas / t+, 
ye a Tas /Toa 
mata f, igh T To4 
SW T - pA Pod ered 
_ arse: mg VTS4 VT 

‘sadiedl 7 Tat ‘FOS OSIGN 
XxNT = BE AE). 

AT = 4. 9006 


BT ’ “Ss 2. + (i= 





_ 
7 
Z 


Py. 


CT = ney 
APs 
my + 
or = waft (a-n) | 
ea 
= 18 = oe — .SCXNTY 
= Sa 
ey = <3 
a a 
aT. = ars - (PRT) | 
“a 
- era 
EFET ‘ [4 ss =a 
GET 
Ter = , — EFFT |i - (Per) | 


Additional Machine Assumptions 


i Fas = Be = Fu 


(no duct loss) 


(no combustor losses) 


Fox = Fes 
The machine characteristics are combined with the 
differential equations to produce a set of system state 


equations. The expanded equations appear in Appendix II. 
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III. THE OPERATIONAL AMPLIFIER 


A. INTRODUCTION 
The fundamentals of the operational amplifier, (OA) 
are briefly presented to establish the necessary analog 


"tools" to proceed to a simulator design. 


B. BASIC OPERATION 


An ideal OP amp can be modeled by the following 


sketch: 


SIGNAL 





with no differential voltage between Rositive a) and 
ReeeCN ) inputs, both emitters are equally biased, and 
collector currents Qu. and Qzare equal. The differential 
voltage output between 1 and 2 is also zero. Driving vr 
more positive with a signal will cause the current through 
Q, to increase. The current through ary must decrease by a 
proportional amount, so that in the differential voltage 
output between 1 and 2, 1 gets driven more negative with 


respect to its initial voltage, and 2 becomes more positive. 





== 


Thus for a small signal applied to one of the input terminals, 
a high voltage gain and current "buffering" takes place as 
a result of the transistor characteristics. 


With the ideal OA, the input current is essentially 


zero, and the voltage gain is assumed infinite. 


c. THE INVERTING OA 


Looking at the following diagram, 





1 


It is assumed that the impedance of the amplifier 


Ze 0 , therefore lez 0. 


Closed Loop Gain, froma summation of currents: 


H 
0 
\l 


C Vo - Va ) /Zz 
: ry ed 
Ta = ONi -— Va) /@i . Tat -Te 


: Ye /\, 
- + Val Zo = Vi fey 2 Sie 


fh 
- : 
2 TN see -/ Cae ve) 


— 4 ‘La 21 /(e ‘22) | 
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If a, the amplifier gain, is large, 

Ne ~~ ~-#2 

ay - ee 
It can therefore be seen that the inverting OA inverts 
and provides a gain proportional to the input and the 


Signal feedback resistances. 


D. INVERTING SUMMER 
f 
By joining multiple input currents to the input side 
of the OA, an inverting summer can be constructed with 


variable gains on each input. 


Vi ZeL 

d Aes No 
V2 ez 

oY 


By current summation, 


No / — Vs. /- ai = Vz. / (at 





aap — 


EB. INTEGRATOR 
Adding a capacitor to the feed back loop of the OA 


adds the following constitutive relation to the feedback 


loop: 








Ta = -Le 
ANo . Eb 
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F. LOGARITHMIC GENERATOR 

Using a diode in the feedback loop enables the 
logarithmic nature of the diode voltage - current relation 
to be included in the loop. The diode characteristic 


can be modeled by the following equation: 


wr 
Z = ew (%)-2 
-——<—> Se ee | =: 
= Lalas — 
Wy Ly 
V = Mo Lu (In (co) 
where: 
Lo = junction saturation current 
Vo = junction saturation voltage 
J 
bel = bias current 
WV = bias voltage 


A typical characteristic appears as: 
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It is important to note that the saturation voltage No is 
directly proportional to the absolute temperature. 


A typical log circuit can be constructed as follows: 


C2 


el Dy, 
: ey 
Vi AKAN A Vo 
Re, 
Sou 
A diode - connected transistor is used to perform 


the diode function, as transistor emitter-base junctions 
have more dynamic current range than diodes in general. 

By the following diagram, an input Vv across Ki causes 

a current flow through the "transdiode". The bias voltage, 


appearing as No varies as the log of “. 
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G. ANTILOG GENERATOR 
In a similar fashion, input voltages or currents may 


be exponentiated in the following circuit: 
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IV. ANALOG COMPUTATION 


A. INTRODUCTION 

To implement the non-linear system equations the 
following functional operations are performed: addition, 
subtraction, integration, multiplication, division, 
exponentiation to a fixed power, and exponentiation to 
a variable power. The first three analog operations have 
been briefly discussed. To implement the higher functions, 


the use of logarithmic circuits must be made. 


Bin THE PRACTICAL LOG CIRCUIT 

There are at least three major considerations in 
the construction of a practical log or antilog circuit: 
cost, complexity, and accuracy. Three different 
architectures are considered to accomplish the log function. 

The first choice is the previously discussed simple 
log/fantilog circuit. Using low quality components, cost 
of the log circuit is under one dollar. Thermal stability 
is considered poor. A test circuit and transdiode 
characteristics are shown in Appendix III. Test data 
indicated that the accuracy and thermal properties of 
this circuit were marginal. 

A second choice is a compensated log circuit as 


diagrammed below, which offers increased thermal stability 
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with the use of a thermistor:14 Nee 
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In the diagram, the collector of 92 is kept at virtual 
ground by the OP amp, and the collector current through 
Oz is equal to a constant reference current provided 
by Veer - It can be Shown that the output voltage is of 
the form NES = ~LOGyzo (zz) where the thermistor 
provides temperature compensation. One of the drawbacks 
to this scheme is that the number of components necessary 
to generate the log function has gone from 3 to at least 
9, with proportional space and cost increases. This 
configuration was not evaluated. 

The third and most expensive choice is to utilize 
a dedicated logarithmic device. The particular unit 
tested was a Texas Instrument log amplifier. This is a 


dual log generator on a single DIP package, at a cost of 
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approximately $2.50. With appropriate circuit manipulation, 
the desired log/anti-log operations can be performed. The 
TI log amplifier was evaluated in depth. Results are 
documented in Appendix III. The log amplifier was the 
foundation of the simulator design. Other dedicated log/ 
antilog modules exist which offer increased dynamic range, 
temperature stability, and accuracy, but alternative 
devices were an order of magnitude more costly than the 


TI log amplifier. 


e. EXPONENTIATION - ALL CASES 

In order to achieve either fixed or variable power 
exponentiation, the use of logarithms is basic. The 
signal flow for a fixed power exponent (power function 


generator or PFG) circuit is diagrammed below: 
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By a judicious selection of values for k, and 2 
current can be multiplied or divided using an operational 


amplifier in the ratio: 
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To be able to exponentiate to a variable power 
requires a cascaded log/antilog arrangement. The use of 
a segment of the state equations is made to demonstrate 


the concept. It is desired to generate the following: 


ZNO = K Cl = ceag=t)) 
where represents the scaling lg several fixed parameters, 
and KW, FARES , [SC and CC are all problem variables. 


The circuit selected must process: 
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As a result of thermal instabilities, the interior section 
of the signal flow was replaced with a four-quadrant 
multiplier. 

The multiplier essentially compresses the separate 


functions into a single chip. By circuit manipulation, 
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a multiplier can also divide, square, and take square 
roots. A range of multiplier devices currently exist, 
with the price of the least expensive about $3. 
Test data on the multiplier used in the simulator 

was collected in Appendix III. To allow as compact a 
circuit as possible, the multiplier was frequently 
utilized in its dividing and exponentiating modes. 
Complete device specifications are included in 


Appendix VI. 
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V. PRELIMINARY DESIGN 


A. INTRODUCTION 

Tina aoe design proceeded in a parallel fashion. 
In an iterative process, estimations were generated for the 
electronic layout, which allowed the cabinet and control 
panel to be sized and designed. 

The simulator is configured to resemble a real-world 
engine control panel. Additional concepts for the per- 
ipheral design which were not incorporated were: 

1. Proportional fuel control schemes with 

variable gains. 

2. Torque-speed loads characteristic of 
pumps, propellers, generators, and other 
real world devices. 

3.  Overspeed alarms and stall audio noise 


circuits. 


B. ELECTRONICS 

imeesoegic Circuit 

The logic processing of the electronic circuitry 
was broken down into 3 separate grid structures, which 
appear as design drawings in Appendix V. Use was made 
of bread-board panels to avoid the inflexibility of other 


techniques. Concurrent with testing various mathematical 
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operations, the circuit design was modified to balance the 
competing parameters of simplicity, size, accuracy, and 
cost. The appropriate equation variables were normalized 
about their design point values, and then scaled to 

provide a circuit signal that was approximately 75% of 

the maximum rates input for the device. The steady 

state design point values for all variables were established 
by a DYSYS program and used for the normalization. 
Additional notes on the wiring techniques and breadboard 


layout are included in Appendix VI. 


Zoe DreYo simulation and Scaling 


A DYSYS simulation of the system equations was run 
as it appears in Appendix IV. Significant comments on 
the program are: 
1. The conditional value of FC was changed 
to a fixed value of 3.0 to facilitate the 
Gireuit logic. 

2. Markunas estimated a parabolic "pump-type" 
load torque, which appears as .0000977(3)**2. 

3. The modeling parameter DT is called TT in 

the program to avoid confusion with the 


DYSYS time-step variables, DT. 
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The scaling procedure started with the following 


listings: 


a. Constants 
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b. Machine Characteristics 
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c. Input Variables at Design Point 
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d. State Variables at Design Point 
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The DYSYS simulation was run to verify Markunas's 
Simulations. Use was made of the program to assist in 
adding voltage and current values to the circuit design, 
and to evaluate the effect of fixing the discontinuous 
constant, FC. FC was conditionally assigned the value 
of 2 or 4 in Markunas work. The steady state effect of 
assigning FC as 3 was investigated and adopted. 

The DYSYS simulation was also run off - design 
to generate a test case for the logic circuit. The 
simulation was run at approximately 50% fuel flow and 
steady-state values were generated for all elements of 


the equations. 





=e 


Some difficulty was encountered in determining the 
correct value for the lumped parameter fluid volume. The 
value which appeared to be consistent with Markunas' data 
was 144 ae An approximation based on real world data 
gave 15 ft" aS a more likely figure. Since the volume 
represents a time constant scaling for the fluid capacitance, 
it has no significance for the DYSYS steady state values. 
Changes in the volume assumed can be readily investigated 
in the simulator dynamic response. 

The scaled voltage values and the system variables 


they represent are drawn on the 3 circuit plans in 


Appendix V. 


OF PERIPHERALS 

The main elements of the remaining design work 
included the chasis, power supply and distribution, 
instrumentation and control hardware, and breadboard 
architecture. The design drawings and sketches included 
in Appendix V illustrate the final peripheral design. 
The photograph below shows the internal layout of the 


Simulator. 
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Appendix VI, specifications, contains additional 
data on the peripheral equipment, its wiring, and operating 


parameters. 
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VI. IMPLEMENTATION AND CHECKOUT 


A. PROCEDURE 

To establish the characteristics of the various 
devices and to investigate the methodology of the design, 
10 tests were conducted. The significant results are 
discussed here in the context of the design process. All 


test data is included in Appendix III. 


iS. LOG GENERATION 

Tests l and 2 determined the voltage - current 
relations for the transdiode - connected OP amp, using 
both 2N2222 and 2N2904 transistors. It became apparent 
that the allowable voltage values of the log circuit 
were so narrow that only current variations could be used 


for the signal, over about a 1 1/2 decade range of values. 


Ca POWER EXPONENTIATION 
Tests 3, 4, 5, and 6 were designed to varify the 
fixed power exponentiation circuitry. A workable 
configuration was wired in test 6, uSing a current - 
dividing resistive network. 
Tests 7, 8, 9, and 10 were conducted to validate 
the variable - power exponentiation procedure. The inability 


of the transdiode - connected OP amp to provide a stable 
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cascaded log circuit led to the investigation of a 
4 quadrant multiplier in Test 8. 

The Exar 2208 multiplier characteristics were 
investigated to generate hookup diagrams, trimming 
values, and allowable voltage values. 

With Test 9, a simple log/multiplier/simple 
antilog circuit was wired and observed. Stability was 
Significantly improved, but was still considered below 
the accuracy level desired. 

In Test 10, a Texas Instrument 441 log amplifier 
was tested to determine its log/antilog characteristics. 
The log-amp was incorporated into the variable expon- 
entiation circuit and the stability evaluated. The 
steady state accuracy of the log amplifier in the 
variable exponentiation circuit was somewhat better than 
the simple log configuration. Since the dedicated 
log amplifier offered better compensation characteristics 
than the simple log configuration, it was substituted 


at all points in the circuit design. 


D. GENERAL CIRCUIT DESIGN 

Concurrent with the test sequence, the overall 
circuit design was changing as new test data was generated. 
The more significant problems which caused iteration of 


the circuit are briefly discussed. One immediate limitation 
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exists uSing the TI log amplifier module. Looking at 

"TL 441 Transfer Characteristics" included in Appendix III, 
the log amplifier only has a useable output range of 0 - 
450 millivolts. This can be doubled by wiring an OP amp 

to the input side of the log amplifier to use the dual 
input capability of the device. Still, the output accuracy 
is limited, since various offsets for multipliers and OP 
amps can easily reach 10 millivolts. 

As a general rule, it was decided to trin the 
circuit at a limited number of points, rather than dealing 
with the offsets required for each component. 

It was intended in the final circuit design to 
operate all components at one dual-polarity supply voltage. 
This choice nominally would be + 8 volts, the maximum 
rating of the most sensitive device, the TI log amp. Test 
10 demonstrated that significant nonlinearities were intro- 
duced into the multiplier characteristics if the supply voltage 
were reduced. Three solutions were postulated: 

1. Run the multipliers at + 15 volts, and the 

OP and log ampsat + 8 volts. 

2. Compensate the multiplier to have linear 

characteristics at + 8 volts supply. 

3. Replace the multiplier modules with log 


amplifiers to process the same function. 
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The first choice was selected in an effort to keep 
the circuit design as simple as possible, in as much as 
the single - voltage supply was more of a packaging 


consideration than an electronic one. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. FEASIBILITY OF THE DESIGN 

From the test data collected, a real time simulation 
can be conducted with the logic circuits wired as designed. 
The techniques developed for the implementation of non- 
linear differential equations should be applicable to a 
wide range of system models. 

Markunas defined his engine in essentially non 
dimensional or pseudo-dimensionless groups. It was not 
possible to investigate the ability of the simulator to 
represent different geometric machines, however, it 1s 
felt that future work will verify the utility of the 


Simulator and the nonlinear model. 


B. DESIGN REFINEMENTS 

The simulator as designed could be significantly 
enhanced with more general torque and fuel control 
schemes. If the simulator could demonstate sufficient 
generality in modeling the behavior of real world gas 
turbines, an autonomous "observer" circuit could be 
configured for a wide range of engines. This observer 
device, using the logic developed in the simulator, could 
be an extremely compact unit. It could provide improved 
control of fuel, speed, temperature, and pressure 


dynamics. Such a circuit could be incorporated in one 





=50— 


or two monolithic LSI chips, to be externally trimmed for 


a particular engine or operating condition. 
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APPENDIX I 


MACHINE CHARACTERISTICS 
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APPENDIX II 


STATE EQUATIONS 
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APPENDIX III 


Test Data 


A. TEST 1 AND 2 


From reference data, the transdiode voltage current 


characteristics was estimated for various values of 


saturation current. The log circuits diagrammed subsequently 


for the 2N2222 and 2N2904 transistors demonstrated low 
base-emitter saturation currents, thus making the exponential 


characteristic curve more abrupt. 
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B. TESTS 3 - 6 

The fixed exponentiation of a variable was evaluated 
in Tests 3 - 6. Tests 3, 4, and 5 were failures as the 
essential concept of current division to achieve the 
exponentiation was not understood. Test 6 demonstrates the 
use of a resistive network which accomplished the current 
manipulation. It is now apparent that as a result of the 
characteristics of the simple log configuration, current 
division can be easily accomplished by assigning the 
appropriate resistor value to the antilog input. This 


is diagrammed: 
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TEST 6 


CIRCUIT DIAGRAMS AND DATA 
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C. TEST 7 


It is desired to validate the circuit segment for the 


(1 — (Pec c-t)) 


iZ2c 


operation 


The operational sequence is as poo ONS 
ANTI - Ce 


ef eH 


Co 


The circuit described in the test data was set up to 
evaluate the feasibility of the operation. 
The desired operating characteristic of the cascaded 


logger - antilogger is picutred as follows: 


ae £2 we K-82 





{D1 oP\. 





agg 


The first log geneator establishes a voltage - 
current relation based on the diode and 20K resistor. 
Consider the second log generator to be superimposed in 
the first with a 3K resistor and a second diode. 

The cascaded antilog - configured amplifiers 
displayed gross thermal sensitivity. The simple OP amp - 
transdiode configuration is unsuitable for this cascaded 


operation without additional temerature compensation. 
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1. TEST 8 
ieee aupose 


Test 8 was conducted to verify the variable exponent 
part of the analog circuitry, using the four quadrant 
multiplier. 


The eguation of interest is written below: 


M-. <= ie (25:15) 
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ii. Methodology 


The functional method used to exponentiate to a 
variable power was to regulate the multiplier output 
current in the ratio called for by the exponentiation. For 


| example: 
; 


d/o te 








= ea 


For multiplier scaling the integral buffer OP amp 
in the 2208 was used for the current division with a 


resistive bridge as follows, 


MULTIGIER + Ruckez OA ASTROS 





such that 1, = i, and the output voltage of the multiplier 
is thus controlling the current input to the antilog 
amplifier transdiode. 


The multiplier output is nominally: 


Meo= Nx Ny 
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It was determined that a resistive bridge must be 


used for the multiplier input from logging OA, sketched 
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as follows: 
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so by modulation of Se input, the desired current ratios 


across multiplier can be achieved. 


iii. Data 
A 5 minute observation of the steady state output 
voltage was made to investigate the circuit thermal 


stability. 
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E. TEST 9 
a Purpose 


To investigate the Texas Instrument TL 441 log 


amplifier. 


1i. Log Configuration 
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The output characteristic VWYY = = Law Ay + Law Az 
verified. From the plot, "TL 441 transfer characteristics, 
the log amplifier has 1 1/2 decades of dynamic input range, 


from 20 - 600 mV. 


iii. Antilog Configuration 


TWO ALTERNATE WIRINGS WEG 
INVESTIGATED ¢ 
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The antilog data is graphed in "TL 441 antilog 
characteristics" for configuration 1. Configuration 2 will 
expand the antilog dynamic range somewhat above the 1 1/2 


decades available for configuration l. 
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F. TEST 10 
fa EUrpose 


To incorporate the 441 log amplifier and the 
quadrant multiplier in the variable exponent circuit. The 
general wiring is sketched below: 


oe 


i rieast 





W/ 


Pieoceaayestate Perrormance 


The circuit demonstrated the ability to track a 


desired variable exponent. It was of primary interest 
to determine if the thermal properties were improved over 


the simple log-fed multiplier circuit. 
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The following data was generated: 


Variable Exponent Test 
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This represents a 50% reduction in the drift over 
the simple log configuration. Since the 441 log amplifier 


is fairly cheap, it is substituted for the simple log 


configuration. 
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APPENDIX V 


DESIGN DRAWINGS AND SKETCHES 
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APPENDIX VI 


SPECIFICATIONS 


ELECTRONIC COMPONENTS 


1. Operational Amplifiers 


Ae 


Fairchild “A 741 


Maximum ratings 
supply voltage 


Power dissipation 
Supply current 


Typical input offset 
voltage 


Fairchild “WA 324 Quad OP amp 





Maximum ratings 
supply voltage 
Single polarity 
Input voltage range 
Power dissipation 


Supply current 


Typical input offset 
voltage 


2. EXAR 2208 Four Quadrant Multiplier 


Maximum ratings 
supply voltage 


Power dissipation 
Typical supply current 


Typical output nonlinearity 


+ 18v 
500 Mw 


2.8 MA 


1 mv 


32vV 
-~0.3v to 32v 
900 MW 


2 MA 


2 mv 


i 18Vv 
750 MW 
5 MA 


5% 
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3. Texas Instruments TL 441 Logarithmic 


Amplifier 

Maximum ratings 

supply voltage + 8v 

Input voltage 6V 

Power dissipation 500 MW 

Typical differential 

output offset voltage 40 Mv 

Typical supply current vec" 18.5 MA 
voc 8.5 MA 


4. Trimming Potentiometers 
All applications: Spectrol type 64 
ceramic potentiometers, 1-100K 

5. Resistors 


All applications: Corning 2% 
1/4 watt resistors 


6. Power Supply 


Analog Devices model 975 dual polarity 
power supply 


Output voltage agen bei, 

Maximum rating 500 MA 

Typical output error 1% 
PERIPHERALS 


1. Breadboard Panels 


Continential Specialties model QT 59 
6.5" panels 





aa 


Volt Meters 


Beede Instrument 3" panel meters, 
model 3-03-8 


0 
250 movement 


20082 internal resistance 


Joystick 


2 - axis movement;two, 100K rotary 
potentiometers on each axis. 
Chassis 


.062" clear anodized aluminum assembled 
with pop-rivet construction. 
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